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ABSTRACT 

We investigate the process of formation and subsequent evolution of promi- 
nence plasma in a filament channel and its overlying arcade. We construct a 
three-dimensional time-dependent model of an intermediate quiescent promi- 
nence suitable to be compared with observations. We combine the magnetic field 
structure with one-dimensional independent simulations of many flux tubes, of 
a three-dimensional sheared double arcade, in which the thermal nonequilibrium 
process governs the plasma evolution. We have found that the condensations 
in the corona can be divided into two populations: threads and blobs. Threads 
are massive condensations that linger in the field line dips. Blobs are ubiquitous 
small condensations that are produced throughout the filament and overlying 
arcade magnetic structure, and rapidly fall to the chromosphere. The threads 
are the principal contributors to the total mass, whereas the blob contribution 
is small. The total prominence mass is in agreement with observations, assum- 
ing reasonable filling factors of order 0.001 and a fixed number of threads. The 
motion of the threads is basically horizontal, while blobs move in all directions 
along the field. The peak velocities for both populations are comparable, but 
there is a weak tendency for the velocity to increase with the inclination, and the 
blobs with nearly vertical motion have the largest velocities. We have generated 
synthetic images of the whole structure in an Ha proxy and in two EUV chan- 
nels of the AIA instrument aboard SDO, thus showing the plasma at cool, warm, 
and hot temperatures. The predicted differential emission measure of our system 
agrees very well with observations in the temperature range log T = 4.6 — 5.7. 
We conclude that the sheared-arcade magnetic structure and plasma behavior 
driven by thermal nonequilibrium fit well the abundant observational evidence 
for typical intermediate prominences. 
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Introduction 



Quiescent solar prominences have been studied for a long time, and significant progress 
in determining their nature has been made. We can say with certainty that they are com- 
posed of cold plasma, and it is generally accepted that the magnetic field supports these dense 
structures in the tenuous corona. However, there are still important unknowns about their 
structure, formation, and evolution. High-resolution observations from space and from the 
ground - e.g., from SOHO, TRACE, SVST, VAULT, Hinode, and SDO - show a complex 
fine structure in space and time. This fine structure co nsists of moving blobs and bun- 



dies of horizontal threads with counterstreaming motions ( 


Enevold 


1976; 


Zirker et al. 


1998 


Lin et al. 


2003; 


Vourlidas & Korendvke 


2003; 


Lin et al. 


2005; 


Bereer et al. 


2008 


). The quies- 



cent prominences are immersed in structures even larger than themselves. Prominences form 
abov e a polarity inversion 



19721 ) in filame nt channels (IMartin 



ine (iBabcock fc Babcock 



magnetic field (iHagyard et al.lll984j : iMoore et al. 



199 



1955 



Gaizauskas et al 



1987 



Smith fc Ramseylll967l : iMcIntosh 



20011). which are regions of sheared 



Venkatakrishnan et al.lll989[ ). Mag- 



netic loop arcades ov erlie filament channels, s 



j anning high a bove them and rooted at both 



sides of the channel (jTandberg-Hansseru Il995l : iMartinI Il998l ). Direct measurements of the 



magnetic fields have shed some light on the rnagnet ic structure of quiescent prominences 
( iLeroyI Il989l : ICasini et al.l l2003l : iKuckein et al.l 120091 ). However, these measurements are 
focused on the cool prominence plasma alone and do not establish definitively the global 
structure of the magnetic field in the coronal portion of the filament channel. 

There are basically two models of the magne tic structure of filam e nt channels: the 
sheare d arcade and the flux rope (see r eview bv iMackay et al.l l2010l ). lAntiochos et al. 
fll994l ). IPeVore fc AntiochosI fcoool ). and lAulanier et al.l f l2002h showed that the sheared 
arcade model exhibits most of the observed properties of the large-scale filament channel 
structure. The generated magnetic field st ructure has d ipped and flattened field hues con- 



sistent with most prominence observations (jMartinlll998|). Manv quiescent filaments are ob- 



served to grow by linkage between smaller segrn ent s (iMalherbe 



1989 



Gaizaus 



<:as et al. 



2001 



Schmieder et al.l l2004at Ivan Ballegooijenl l2004j ) . iDeVore et al.l (|2005[ ) and lAulanier et al. 

studied the interaction of two neighboring sheared bipoles with different combina- 
tions of polarity and chirality. They found that merging two adjacent sheared arcades with 
aligned axial field formed a long, stable, strongly sheared prominence. Assuming that the cool 
plasma resides in magnetic dips, the resulting structure reproduces the glob al appearance 
of a typical low-latitude intermediate quiescent filament (jMackay et al.ll2010l ). The present 
work goes beyond this idealized assumption to determine the origin and time-dependent 
evolution of the cool prominence plasma. 



It is generally believed that the origin of the prominence plasma is the chromosphere 
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( iPikernerlllQTll : ISaito fc Tandberg-Hanssenlll973l : IZirker et al. I ll99i ). 3ut the mechanism that 



puts the chromospheric material into the corona is under debate. O ne promising candidate. 



the thermal nonequilibriuin process, has been studied extensively ( Antiochos &: Klimchuk 



199ll : lAntiochos et al.lll999L I2OOOI : iKarpen et al.ll200ll . I2003L l2005l . I2OO6I : iKarpen fc Antiochos 



20081 ). and yields quantitative predictions of prominence plasma directly comparable with 



observations. In this model, localized heating above the flux tube footpoints produces evap- 
oration of the chromospheric plasma, which condenses in the coronal part of the tube. 
The localized heating near the chromosphere is c onsistent with several coronal-loop studies 



( lAschwanden et al.ll200ll : ISchmieder et al.l l2004b( ) . To produce this instability the heating 



must be concentrated in a region of approximately one eighth of the tube length or less 



fISerio et al.lll98ll : iMok et al.lll990l : lAntiochos fc KlimchukI Il99ll : lOahlburg et al.lll998r ). In 



our earlier investigations of the thermal nonequilibrium process, we found that large conden- 
sations form along long, low-lying magnetic field lines. In most cases the plasma undergoes 
repeated dynamic cycles of plasma condensation, streaming along the tube, and falling to the 



Schrijver 


2001; 


Miiller et al. 


2003. 


2004, 


2005) 



condensations form readily but the subsequent evolution depends on the tube slope. For 
slopes steeper than a critical value determined by the heating asymmetry and gravitational 
scale height, the condensations fall to the bottom of the dip and remain there, accreting 
mass continuously. 

In the present work, we built a comprehensive model of a multi-threaded p rominence 
combi ning the results of the three-dimensional (3D) sheared double-arcade model of iDeVore et al, 
(120051 ) for the magnetic structure with the thermal nonequilibrium process for the plasma 
evolution determined by many independent one-dimensional simulations. The resulting 
three-dimensional model allowed us to study formation and evolution of the plasma in a 
multi-stranded prominence. In §5] we describe the theoretical model used in the current 
work. The detailed results of the investigation are presented in §3J The global appearance 
of the prominence in different spectral channels and the predicted differential emission mea- 
sure are shown in §H Finally, in ||5] the results of this investigation are summarized and 
conclusions are drawn. In particular, we conclude that our three-dimensional model fits well 
the observational evidence. 



2. Theoretical model 



As in our earlier studies, we a ssume that the magnetic structure is static on the timescale 
of the plasma evolution (see e.g., lMartinlll998[ ) and that the gas pressure of the plasma is 
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negligible in contrast with the magnetic pressure. Additionally, in the solar corona the 
conductive heat flux is mainly constrained to follow the field lines. With these consider- 
ations each flux tube is a thermally isolated system where the plasma motions are along 
the magnetic field. The equations describing the plasma evolution are reduced to one- 
dimensional (I D) set of hydro dyn amic equations for mass, momentum and energy conserva- 



tion (see. e.g..lKarpen et al. 



(IKlimchuk fc Carrillll2001f l. 



20051) includi ng the Raymond- Klimchuk radiative loss function 



Spitzerl ( 1l962l ) thermal conductivity, and a prescribed volumet- 



ric heating. Hence the time-dependent plasma evolution along each static field line can be 
determined by solving numerically those equations. 



We use the magnetic field structure given by the sheared double-arcade model of lDeVore et al. 



(I2OO5I ) in Cartesian geometry. This magnetic field configuration is the result of shearing two 
identical dipoles with the same polarity, the same chirality, and with their axes aligned. The 
shear is produced by imposing differential footpoint motions over the photospheric plasma 
near the polarity inversion line (PIL). As a consequence, an inner bundle of sheared field lines 
is produced (see Fig. [1]). These are long, low-lying field lines nearly parallel to the PIL, and 
many of them have dips. An outer bundle of un-sheared field lines comprises the overlying 
arcade with arched field lines perpendicular (or slightly skewed from perpendicular) to the 
PIL. In this model, the filament channel (FC) magnetic configuration is a natural conse- 
quence of the magnetic shear in a 3D structure. However, the exact origin of the magnetic 
shear is not important in our model. The shear could be the result of the partial emergence 
of twisted flux tubes from below, or direct subsurface motions near the PIL. 

We performed the numerical simulations of the one-dimensional hydrodynamic equations 
with our code ARGOS (Ad aptively Refined GOdunov Solver) with adaptive mesh refinement 



see 



Antiochos et al.lll999l for details). We have selected a discrete set of 125 non- uniform 



flux tubes of the 3D sheared double- arc ade magnetic s tructu re shown in Figured! We have 
rescaled the computational domain of iDeVore et al.l ( 120051 ) to a box of 800 Mm parallel 
to the PIL, 200 Mm perpendicular to this PIL in the base plane, and 200 Mm in the 
vertical direction. The footpoints of the field lines were selected by imposing two uniform 
rectangular grids in the photosphere. The first grid is denser and covers the sheared region 
of the photosphere, being 40 Mm wide centered at the PIL and 600 Mm along the PIL. The 
second grid covers the outer region with no shear, and it is located only on one side of the 
PIL due to the symmetry of the magnetic field with respect to the PIL. This grid extends 
perpendicular to the PIL from the end of the first grid to 55 Mm away. Along the PIL both 
grids have the same extension. The geometry of the 3D flux tubes is incorporated into our 
ID simulations in two ways. The field line shape determines the projection of the gravity 
along the tube. In addition the non-uniform cross-sectional area A is incorporated using 
conservation of magnetic flux along the tube. The position s along the tube is described by 
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Fig. 1. — Magnetic field structure of tlie slieared double arcade in Cartesian geometry. White 
lines are the magnetic field lines we have used in this work. In (a) the magnetic structure is 
viewed from above. At the bottom of the corona the normal magnetic field is plotted (color 
scale). The two dipoles are clearly seen in this view. In (b) the side-view orientation is 
shown. Finally, in (c) the end-view orientation is plotted. 

the distance from one footpoint at the base of the corona. With this definition the positions 
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Fig. 2. — Examples of the three types of flux tubes in our prominence model. The horizontal 
d axis connects the two tube feet. Flux tubes with a deep dip (solid lines) are sheared lines, 
and form the filament channel core. Shallow-dip tubes (dotted lines) are also sheared lines, 
and cross the polarity inversion line above the deep-dip lines. Finally, the loop-like tubes 
(dot-dashed lines) have no dip, and form the overlying arcade structure. 

of the two footpoints at the base of the corona are s = and s = L, where L is the length 
of the coronal part of the tube. Setting the area at s = to be Aq, we can determine the 
area along the entire tube using A{s)B{s) = AqBq where Bq is the magnetic field strength 
of the sheared double arcade at s = 0. 

The geometry of the flux tubes in the whole magnetic structure can be classified roughly 
into three groups: tubes with deep dips, tubes with shallow dips, and tubes with no dips 
(loop-like tubes). In Figure |2] two examples of each group are plotted. The first group has 
deep dips; these are sheared, long, and low-lying lines with inverse polarity at their central 
parts (tubes 1 and 2 in Fig. [2]). The magnetic field and gravity jointly provide a stable 
equilibrium region in the dipped portion of the tube because the projected gravity forces the 
plasma to move to the center of the dip. The second kind of tube has a plateau shape with 
a central part almost horizontal, as we see in the examples 3 and 4 of Figure [21 providing 
a region in the tube where the projected gravity is small or zero. These two groups of field 
lines belong to the sheared magnetic field, and they form the FC structure (Fig. [1]). The 
third kind of tube is in the unsheared or weakly sheared arcade overlying the FC. These 
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tubes are loop-shaped as exemplified by tubes 5 and 6 in Figure [2l Most of the loops cross 
the PIL perpendicularly, but the short tubes near the FC structure are skewed from the 
perpendicular direction. Most loop-like tubes have a large projected gravity except in the 
tube apices. There is a small group of loop-like tubes at both ends of the FC with large 
inclinations (see Figs. [1^ and[T]D), with small projected gravity along the tube. Tube 5 is a 
short loop that wraps around the inner core of the FC structure and is slightly skewed from 
perpendicular to the PIL. Tube 6 is a very long loop reaching very high altitudes of around 
190 Mm. Most of the selected 125 tubes belong to the FC, and the others belong to the 
overlying arcade. 

In Figure [3] the histogram of the lengths of the flux tubes is plotted. The distribution 
has a Gaussian shape except for the contribution of very long tubes. The tube lengths are 
distributed between 84 Mm and 479 Mm with a mean value of 218 Mm. We compute also 
the contribution of the three kinds of flux tubes to the global distribution of tubes. We 
see that the right wing of the distribution (i.e., L > 218 Mm) is given mainly by the tubes 
of the first kind with deep dips; the main contribution to the left wing of the distribution 
comes from both shallow-dip and loop flux tubes. Thus, on average the tubes with deep 
dips are longer than the tubes with shallow dips. This is because lines in the first group 
are reconnected during the process of formation of the FC by interactions between the two 



sheared bipoles (iDeVore et al.ll2005l ). We have selected only a small number of tubes of the 
overlying arcade, as can be seen in the histogram. This is because we are interested mainly 
in the tubes of the FC core, but we also included a small sample of other tubes to take into 
account the evolution of the plasma outside the FC. The area expansion factors, A{s)/Aq, 
of the tubes of the first and second kinds in the corona are similar, with a mean value of 
approximately 3 in the central part of the tubes. In contrast, the area expansion of tubes of 
the overlying arcade has a mean value of 30 in the loop apices. 

As in our previous studies, the volumetric heating has two components: a spatially 
uniform background heating, Hq, plus a localized component at each footpoint in the chro- 
mosphere. Hi and H2, that falls off exponentially in the corona. Thus, 

H{s) = Ho + Hi e^'/^ + H2 e(^-^)/\ (1) 

where A is the heating deposition scale that we have set to A = 10 Mm, and the background 
heating is set to Ho = 1.5 x 10~^ erg cm~^ s~^. The tube feet are heated asymmetrically, so Hi 
and H2 are different. We have set the larger value of the volumetric heating as max(ifi, H2) = 
10^^ erg cm~^ s~^, and the heating asymmetry as mm{Hi, H2) /in.eix{Hi, H2) = 0.75. In 
our 125 independent simulations we have distributed randomly the position of the larger 
localized heating, so some tubes are more heated at s = and others at s = L. A 60— Mm 
chromospheric region has been added to each flux tube footpoint in order to have a mass 
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Fig. 3. — Histogram of the lengths of the selected flux tubes of the magnetic structure. The 
full set of 125 (black dashed line) has a distribution similar to a Gaussian. The minimum 
length is 84 Mm and the maximum is 479 Mm, with a mean value of 218 Mm. Tubes of 
the first kind (blue line) with deep dips have a Gaussian distribution with a mean length 
of 255 Mm. The second kind (green line), tubes with shallow dips, also have a Gaussian 
distribution but centered at 156 Mm. The third kind (red line), the overlying arcade loops, 
are tubes with a mean length of 192 Mm. 



reservoir and sink for the coronal plasma evolution. In all flux tubes, an equilibrium situation 
was established first with only background heating Hq turned on during the first 10^ s. This 
small initial uniform heating p opulates the flux tub es with density, temperature, and pressure 
dictated by the scaling laws of lRosner et al.l ( 119781 ). After this 10^ s, the localized heating was 
ramped up over 1000 s and sustained at that level thereafter for an additi onal 10^ s. Further 



details of the numerical technique can be found in lAntiochos et al.l (Il999[ ) and iKarpen et al. 

fcoo5h . 



3. Results 



The evolution of the plasma o n each flux tube de pends on the tube geometry and the 
asymmetry of the heating function (IKarpen et al.ll2003l ). As we have seen in §2] the geometry 
of the fleld lines can be classifled roughly into three groups. In Figure |21 two representative 
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tubes of each group have been plotted. In order to understand the process of formation 
of the whole prominence, we first study the evolution of the plasma on the selected set of 
six field lines. In Figure IH we have plotted the mass of the condensations as a function of 




Fig. 4. — Temporal evolution of the condensation masses per unit area for the six represen- 
tative tubes, where the masses are divided by Aq. (a) Condensations in the two tubes with 
deep dips, (b) Condensations formed in tubes with shallow dips, (c) Condensations in the 
two overlying arcade tubes. 
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time for the six simulation examples. Figure shows the evolution of cool condensations 
in tubes 1 and 2 with deep dips. The condensation in tube 1 forms at t = 3 hrs after the 
onset of the localized heating, after which the condensation mass grows linearly at a rate of 
6.6 X 10~^ g cm~^ hr~^ for the rest of the simulation. The condensation forms at the dipped 
portion of the field line and remains there while it accretes mass. Each condensation of tube 
2 shows completely different behavior, formin g close to a footpo i nt and falling rapi dly to 



it. This process is repeated in a limit cycle ( iKarpen et al.l l2003l : iMiiller et al.l l2003l ) with 



a period of ~ 3 hrs. In this cycle a pair of consecutive condensations forms, moves along 
the tube, and finally falls to the chromosphere. The difference between tubes 1 and 2 in 
terms of the plasma evolution is dictated by t he heating asymmetry relative to the flux tube 



geometry (see, e.g., iKarpen et al.ll2003l . |2005| ). In both cases the heating is stronger in the 



righthand footpoints of Figure [21 In tube 1 the steep slope between d = 40 Mm and 110 
Mm keeps the condensation in the dip. However, in tube 2 the slope is not steep enough to 
retain the condensation, so it falls to the lefthand chromosphere. 

The evolution of the mass in shallow-dip tubes, simulations 3 and 4, is shown in Figure 
I3)d. In these two simulation cases there are no stable condensations as in simulation 1; only 
small and short-lived condensations are produced cyclically (note the change in scale between 
panels of Fig. H]). In simulation 3, the sequence of condensing and falling to the chromosphere 
is repeated every 2.5 hrs, and each condensation lives for approximately 0.5 hrs. However, 
simulation 4 produces more massive condensations roughly every 4 hrs, each of which has a 
lifetime of 3 hrs. Tube 3 has a length of 176 Mm, and tube 4 has a length of 240 Mm. This 
indicates that the longer the loop, the larger the period of the limit cycle. 

Similarly, the evolution of the cool mass of the two loop-like tubes is plotted in Figure 
Ht. Small condensations form in both loops of the overarching arcade. These two tubes (5 
and 6) have lengths of 121 Mm and 479 Mm respectively, with tube 6 being the longest tube 
considered in this work. In simulation 5, we see that each condensation forms quickly, lives 
for ~ 0.25 hrs, then falls rapidly to the chromosphere. This process is repeated periodically 
every 3 hrs. The temporal evolution of the cool mass in loop 6 is more complex, with several 
peaks reflecting a cluster of short-lived condensations that forms and disappears within an 
interval of ~ 6 hours. Only one entire condensation cycle occurs within the simulation 
time, but another condensation appears at the end of the run. Thus, the period of the 
evaporation/condensation cycle is around 18 hrs. 
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Fig. 5. — Plot of the maximum mass of the condensation as a function of tube length. Dia- 
mond symbols correspond to condensations of the Population T (see text) with masses larger 
than the separator mass (dotted line). Asterisks correspond to condensations of Population 
B with masses smaller than the separator. The 6 representative tubes are numbered and 
surrounded with black circles (see text). 



3.1. Condensation masses and lengths 

Figure [5] shows the maximum mass of the condensation in each flux tube. This quantity 
is useful for determining the ability of the tube to produce and accumulate cool mass accord- 
ing to its geometry. When the condensation mass grows monotonically, as in simulation 1, 
this maximum mass coincides with the value at the end of the simulation. When consec- 
utive evaporation/condensation cycles are produced, as in simulations 2 to 6, we average the 
maximum mass at each cycle as Mmax = Yli fnf^^^/N, where m™*^^ is the maximum mass at 
i-th cycle. From this figure we see that the condensations can be divided mainly into two 
separate groups: condensations of relatively large masses (diamonds), which we denote Pop- 
ulation T (for reasons to be explained shortly), and condensations of relatively small masses 
(asterisks), which we denote Population B. It is important to note that all 125 simulated 
flux tubes studied here produce condensations. In the first group, condensations form and 
then grow monotonically until the end of the simulation, as for simulation 1 (see Fig. [2]). 
In the second group, each condensation forms, moves along the tube, and finally falls to the 
chromosphere, as in tubes 2 to 6. We have found that approximately 7 x 10~^ g cm~^ is the 
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separator between the two condensation populations (dotted line). Condensations just above 
this limit initially have cycles of evaporation/condensation and fall to the chromosphere, but 
eventually the condensation mass grows monotonically. Condensations just below the limit 
start to grow monotonically but at some time the entire condensation falls to the chro- 
mosphere, and several cycles are produced in the simulation. All flux tubes that produce 
condensations of the first group are tubes with deep dips, but not all tubes with deep dips 
produce condensations of the first group, as discussed above for tube 2. Small condensations 
of the second group are produced in all flux tube geometries, as the representative examples 
2 to 6 indicate. There is no clear dependence of the condensation mass on the tube length. 
Condensations of the first group are produced in a range of tube lengths between 160 and 
340 Mm. This reflects the range of lengths of the deeply dipped field lines (see Fig. [3]). In 
Figure E] we see a few examples of the second group of condensations for tube lengths larger 
than 350 Mm. These condensations are generated in the overlying arcade, and as mentioned 
in §2] we have only a small sample of those field lines. However, for our small sampling of 
the overlying arcade, the maximum condensation mass does not depend on tube length. In 
summary, all three kinds of tube geometries produce condensations of the Population B, 
with small masses, but a subset of deeply dipped flux tubes yields the massive condensations 
of the Population T. 

In Figure |6l the mass of the hot plasma and the mass of the condensations in all tubes 
are plotted. The mass of the hot plasma (dotted line) grows steadily from a small value 
(25 X 10~^ g cm~^) given by the initial, uniformly heated, equilibrium. In this interval 
the initially underdense tubes are being filled with hot plasma as a consequence of the 
chromospheric evaporation at the tube feet. The growth of the hot plasma slows down at 
about t = 2.0 hrs, and finally stops at t = 2.5 hrs (see inset plot of Fig. [6]). In the interval 
2.0-2.5 hrs, the chromospheric evaporation has increased the coronal density of most of the 
flux tubes and also their radiative losses, so the tubes start to cool down. At t = 2.0 hrs, the 
first condensations appear and the mass of the cool plasma starts to grow. It is interesting 
to note that in the interval 2.5-4.0 hrs the coronal mass is reduced by 30% of its maximum 
value. This partial evacuation of the corona occurs because the cool condensations accrete 
mass faster than the process of evaporation from the footpoints, due to a pressure deficit tha t 
pulls additional material into the condensations from the corona (see iKarpen et al.ll2003l ). 
The growth rate of the condensation mass is higher during this interval than in the rest of the 
simulation. After this time {t = 4.0 hrs) the hot mass reaches a more or less constant value 
of 158 X 10~^ g cm~^, on average, for the rest of the simulation. Similarly, the condensation 
mass reaches a steady growth rate of 39 x 10~^ g cm~^ hr~^. In this steady state, the rate at 
which mass condenses is balanced on average by the rate of evaporation from the footpoints 
of all tubes, so the mass of the hot coronal plasma remains more or less constant. 



- 13 - 




Time (hrs) 



Fig. 6. — Temporal evolution of the mass of the plasma in the filament channel. The mass 
of the condensations (solid line) is separated according to the contributions of Population T 
(dashed line) and Population B (dot-dashed line) condensations. The mass of the coronal 
plasma (dotted line) grows from an initial equilibrium value beginning at the localized heating 
onset (t = hrs). The inset plot shows details of the early evolution. 



We have separated the contributions of the Populations T (dashed line) and B (dot- 
dashed line) to the total cool mass. Figure |6] shows that there is only a slight difference 
between the condensation mass of the entire system and the mass in tubes of Population T. 
This indi cates that themain contribution to prominence mass comes from the Population 
T tubes. iLin et al.l (|2005[ ) estimated the observed width of the prominence threads to be 
200 km. Assuming this value for the average diameter at the deep-dip tube centers, we 
can estimate the total mass of our model prominence t o be 1 x 10^^ g. This agrees with the 
observed values ranging from 1 x 10^'^ g to 2 x 10^^ g (see lLabrosse et al.ll2010l . for an extended 
review and references therein) . The mass estimate depends strongly on the volumetric filling 
factor of the cool plasma within the surrounding hot corona, which can be estimated by 
dividing the coronal volume occupied by cool threads by the total volume in which those 
threads are sparsely distributed. This total volume was determined as described in §1] and 
is ~ 46 X 10^ Mm^. The 47 Population-T threads collectively occupy ~ 54 Mm^, assuming 
cylindrical threads of 200 km in diameter and the average thread length at the end of the run 
derived from Figure [H Therefore the filling factor of our model is ~ 0.001, which matches 
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the lower limit of observed values (ILabrosse et al.ll2010l ). To obtain a larger filling factor, as 
reported for some observations, either the threads are thicker than assumed (which would 
contradict the highest resolution data) or the prominence contains appreciably more threads 
than were selected for our modeling investigation. 
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Fig. 7. — The maximum mass of every individual condensation (see Fig. |5] for symbol 
definitions) versus its length. 

Figure [7] shows the maximum condensation mass as a function of its maximum length. 
We see that condensations of the Population T are long, whereas condensations of the Pop- 
ulation B are mainly short. This result is expected because, when a condensation grows by 
accreting mass, it also grows in size. For the Population T there is no clear relation between 
the mass and length of the condensations. Lengths are in the range 26.4 to 56.6 Mm, with 
a r nean value of 37 Mm. These large condensations form threads similar to those reported 
by iLin et al.l (|2005[ ): the computed lengths are in agreement with the observed lengths. We 
hereafter call the condensations of Population T "threads" . Condensations of Population B 
are smaller, with a mean maximum length of 3.3 Mm. Hereafter, these small features are 
called "blobs" . The advent of high- resolution high-cadence ob servations has revealed that 
such small, dynamic features exist throughout the corona (e-K - - ISchriiverl l2001t iMiiller et al. 
2OO3I : be Groof etldlbooi iMiiller et al.ll2004i Ide Groof et al.lboosh . 



In Figure El the mean lengths of the two populations of condensations are plotted. The 
size of the threads (solid line) grows rapidly between t = 2.0 and t = 3.5 hrs, when the mean 
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size reaches 8 Mm. This period of time coincides with the partial evacuation of the coronal 
mass discussed in §3.11 Note also the sudden increase of the mean size between t = 3.0 hrs 
and t = 3.3 hrs; in this 18— min interval the mean length grows from 2 to 8 Mm. Later on, 
the length remains more or less constant at about 8 Mm, until t ~ 5 hrs when the growth 
resumes up to the end of the simulation. In contrast to the total mass (Fig. [6]), the length 
does not grow linearly, and its growth rate decreases with time. Two factors contribute to 
the diminution of the growth rate, both of which are related to the geometry of the deep- 
dip flux tubes. First, after the formation of the thread the newly accreted plasma occupies 
regions of the tube with larger projected gravity than in the middle of the thread due to 
the concave-up shape of the dip. Then, accreted plasma significantly compresses the already 
condensed plasma because the pressure scale height in the condensation is small (~ 1 Mm). 
A similar result was reported by lKarpen et al.l (120061 ) . who found that only nearly horizontal 
tubes can support extended condensations. Second, the newly accreted material fills regions 
of the tube that are wider than the central portion of the dip, due to the area expansion 
intrinsic to the sheared arcade mechanism for dip formation. In contrast to the threads, 
blobs rapidly reach a constant length of about 1 Mm. These small condensations can form 
anywhere in the tube and are short-lived (see Fig. Hj). This continual process of condensation 
and falling to the chromosphere keeps the temporally mean length more or less constant. 
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Fig. 8. — The mean length of the Population T (solid line), and Population B (dot-dashed 
line) condensations function of time. 
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3.2. Condensation velocities 
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Fig. 9. — Absolute values of the maximum bulk speeds of the condensations of both popu- 
lations as a function of the angle between the bulk velocity vector and the vertical z-axis, 6. 
The symbols are defined as in Fig. [5l The angle ranges from 6 = tt/2 (horizontal motion) 
to ^ = TT (vertical downward motion). 

After formation a condensation moves along its field line in a manner determined by 
the local gravity along the tube and the localized-heating imbalance between the tube foot- 
points. Condensations are extended regions of cool and dense material in which each point 
moves with its own velocity. These motions can be decomposed into a bulk motion of the 
condensation plus local expansions or contractions. Here we investigate the dynamics of 
these individual condensations considering the bulk motion as the mass-averaged velocity 
in the condensation, i.e., the motion of the center of mass. In Figure [9|, the absolute value 
of the maximum bulk velocity of the condensations is plotted as a function of the angle 
between the bulk velocity vector and the vertical direction. Two kinds of motions are asso- 
ciated with the two condensation populations. Thread maximum velocities are concentrated 
in an angle range 6 = 7r/2-2.0 that is almost horizontal. Figure [9] reflects the fact that 
motion of the threads is in the concave part of the tube, and the motion in the dip is mainly 
horizontal with a small vertical projection, resembling a swinging motion. The direction 
of the velocities is not included in the figure, but in the movies associated with Figures [10] 
and [H] we see horizontal counterstreaming in adjacent threads. These condensations form 



-17- 



in the dippe d part of the tube, and fall to the bottom of the dip with a damped oscilla- 
tory motion (lAntiochos et al.ll2000l : iKarpen et al.ll2003l ). The maximum velocity is reached 
beyond the bottom of the dip; this is why the angles are not exactly 7r/2 in Figure O In 
contrast, blob velocities are in the range 6' = 7r/2 — 2.9, indicating that their motions cover 
the full span from nearly horizontal to nearly vertical. This result is expected because blobs 
are produced on flux tubes of all three kinds described in §2l with geometries ranging from 
mostly horizontal tubes to mostly vertical tubes. The motion of the blobs on tubes with 
either deep dips or shallow dips (simulations 2 to 4) is relatively horizontal, while that of 
blobs formed in the overlying arcade (simulations 5 and 6) is nearly vertical. Most of the 
blobs slow down as they approach the tube feet, and in some cases several bounces occur 
before they fall finally to the chromosphere (see movies associated with the Ha prox y of Figs. 



[TUl an d [TT]) . In general the maximum velocities are not reached at the tube feet. ISchrijver 
reported similar motions in the coronal rain observed by TRACE, and suggested that 
the deceleration is the consequence of the p ressure gradients of the underlying plasma close 
to the tube footpoints. iMiiller et al.l ( 20031 ) simulated this phenomenon and confirmed that 
blobs increase the pressure of the plasma below as they approach the chromosphere. Threads 
have a maximum velocity of 83.5 km s~^ and a minimum of 7.8 km s^^, with a mean value of 
34.2 km s^^. Blobs have a maximum velocity of 92.4 km s~^ and a minimum of 16.6 km s~^, 
with a mean of 51.0 km s~^. There is no clear dependence of the bulk velocity on the angle 
in the thread population. However, the blobs show a small tendency for the maximum bulk 
velocity to increase with the angle. This indicates that the downward velocities are larger in 
vertical tubes, but these velocities are far from the free-fall speeds. As we see in Figures [TU] 
and [m the blobs form mai nly below 30 Mm in height. At this height t he free-fall velocity is 
approximately 130 km s~^ (IMackay fc Galsgaardll200ll : ISchrijveii 120011 ). indicating that the 
condensations are slowed down by pressure forces impeding the free fall. 



4. Observational diagnostics 

As in our earlier works on thermal nonequilibrium, we illustrate the dynamic behavior of 
the entire prominence model by computing synthetic images in selected representative EUV 
passbands. Unlike our prior studies of individual tubes seen from a single point of view, 
these images reflect the 3D structure of the prominence model with orientation and line-of- 
sight integration effects included. We have assumed that the emission in the EUV channels 
is optically thin and computed the images using the temperature response function of the 
Atmospheric Imaging Assembly (AIA) instrument aboard the Solar Dynamics Observatory 
(SDO) satellite. We show the AIA band passes 17lA, and 21 lA. Several spectral lines 
contribute to each EUV passband, associated with different ions and ionization degrees. The 
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peak temperature of the 17lA passband is log T = 5.8, and the main contribution comes from 
Fe IX. The 21lA band has a peak temperature of log T = 6.3, and the main contribution is 
from Fe XIV. In addition, a p roxy for the Ha emission is computed as in our earlier papers 



(see, e.g.. 



Karpen et al.ll200l[ ). setting the intensity constant whenever T < 35,000 K and 
zero at higher temperatures. We assume that the flux tubes are unresolved by the telescope, 
and use the AIA instrument's point spread function half-width of 0.6". 

We have constructed off-limb images in two perpendicular orientations of the system 
relative to the line of sight (LOS). In the end-view orientation the LOS is parallel to the PIL 
or filament-channel axis. In the side view the relative orientation of the LOS is perpendicular 
to the PIL. Figures [TOl and [TT] show the temporal evolution in all of the considered channels 
in the end and side views, respectively. In both views different types of structures are visible 
in different emission lines. In the Ha channel we see the temporal evolution of the cool 
condensations, and distinguish between the two condensation populations. As we described 
in §3.11 the threads are the largest cool condensations, containing most of the prominence 
mass. From the Ha images we identify two separated prominence segments in the sheared 
double arcade. Both have a trapezoidal shape with a base of approximately 20 Mm width 
(end view), 115 Mm length (side view), and 20 Mm height. Threads form in an off-center 
position in the dipped portion of the flux tubes, and oscillate around the bottom of the 
dips. In Figure [TTb the recently formed threads are shown. The mean length of the threads 
is approximately 10 Mm (see Fig. |8]), and they move mainly horizontally with a large 
amplitude. In the n aovies accompan y ing t he figures we clearly see the opposite motions 



on adjacent threads. lAntiochos et al.l ( 119991 ) proposed this mechanism of condensation and 



oscillation in the d ips of the flux tubes to explain the counterstreaming motions reported by 



Zirker et al.l ( 1l998l ). Here we clearly see how a bundle of threads manifests counterstreaming 
motions. These motions are also present in the end-view orientation (Fig. ITOb). However, 
the direction of the motion is slightly different with respect to the LOS and the projected 
amplitude of motion is small. In Figure [TTb the threads are larger than in the previous frame 
[TTh . with a mean length of 35 Mm. We clearly see the concave-up shape of the threads in Ha 
(Figures [TOb and [TTb) lying in the bottom of the dips. The amplitude of the oscillations is 
now much s maller than in th e initia l stages of the oscillation, because of the strong damping 



reported by lAntiochos et al.l (120001 ). 



Condensations of the population B (blobs) are small features seen in the Ha synthetic 
images in all regions of the filament channel and the overlying arcade. These ubiquitous 
features form mainly at heights below 50 Mm, although some blobs form at greater heights 
up to 80 Mm. Many blobs condense and move along the field lines between the two thread 
pyramids, producing a counterstreaming flow of blobs in that region (see the accompanying 
movie of Figs. [TTk and [TTb). Blobs forming in the overlying arcade rapidly fall to the 
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Fig. 10. — Time evolution of the whole system in end-view with different filters. In (a) 
and (b) the Ha proxy is shown aX t = 5.56 hrs and t = 23.81 hrs, respectively. Similarly, 
(c) and (d) show these two stages of the evolution in the 17lA passband of the AIA/SDO 
instrument. In (e) and (f) the temporal evolution in the 211A emission intensity is plotted. 
(An animation of this figure is available in the online journal.) 
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chromosphere with an accelerated motion following the magnetic field. As we pointed out in 
§3.2[ these motions are not free-fall movements. When blobs approach the tube footpoint, 
their movement decelerates before they enter into the chromosphere. 

Figures [TOb andfTOH show synthetic images in the 171A passband in the end-view orien- 
tation. In this channel the emission comes from relatively warm plasma at log T = 5.8. In 
these figures, we clearly see strong emission from the core of the filament channel, where the 
prominence resides. We say that the emission from a region or a feature is bright or strong 
when the intensity is larger than three times the standard deviation of the intensities of all 
pixels of the entire image. The source of this emission is the warm plasma in the prominence- 
corona transition region (PCTR) at both ends of the individual threads and blobs (see the 
corresponding Ha images). In this orientation the maximum intensity per image averaged 
over the simulation time is 637 DNpix^^s^^. The emission in this wavelength is strong in 
the core because the threads are concentrated in this region and the bright features are more 
or less aligned with the LOS. 

In the side views (Figs. [TTk and [TTb) the emission from the filament core is not as 
strong, with a temporal average of 374 DNpix~^s~^, because the bright features are less 
concentrated and not aligned with the LOS. However, from this point of view we clearly 
see the bright structures resembling two opposing tadpoles (or comets with tails) on both 
ends of the threads. Between each pair of tadpoles is an emission gap associated with 
the cool material of the condensation. These tadpole features are the consequence of the 
temperature and density distributions at the ends of the condensations. The temperature 
rapidly rises from several thousand K near the condensation to several million K in the hot 
parts of the tube, while the density inversely declines. In an intermediate region between 
the condensation and the corona, the temperature and density conditions are appropriate 
to produce the bright emission. The gradients of the temperature and density steepen close 
to the condensation and fiatten farther away, explaining the tadpole shape of the emission. 
These bright features move together with the condensations; therefore, as shown in the 
accompanying movies, counterstreaming motions are also evident in the 17lA passband. 

The end-view images reveal that the emission from the outer region of the prominence 
is fainter than from the filament core. In this region, the bright features are associated with 
Population B blobs. Since blobs are small condensations, the emission gap between their 
pairs of tadpoles is smaller than in case of the threads. The emission depends on the motion 
of the blobs: the leading tadpole is brighter than the trailing tadpole, relative to the direction 
of the blob motion, and in some cases the trailing tadpole is very faint. The animations also 
reveal long, transient brightenings, which in some cases illuminate most of the tube. These 
very long, bright features are consequences of the evaporation/condensation cycles. Before 
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Fig. 11. — Same as Fig. [TO] showing the temporal evolution in emission intensities in the 
three filters, but from the side-view orientation. (An animation of this figure is available in 
the online journal.) 
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a condensation is produced, a large portion of the material of the flux tube cools down as 
the density increases. When the temperature of this denser plasma is approximately the 
peak temperature of the filter, a strong brightening is produced in that passband. The 
plasma continues to cool down, and finally the condensation is produced. At this moment 
the long, bright feature disappears and the tadpoles appear in 171 A. All of the flux tubes in 
the filament channel produce such long features. However, the tubes with the Population-T 
condensations produce this transient brightening only in the initial stages of the simulation, 
just before thread formation (around t = 2 hrs). Thus, we expect large-scale transient 
brightenings to occur in the 171 A passband, both in the core of the FC and in the overlying 
arcade. 

Figures [TUfe and dUf illustrate the temporal evolution of the system in the end-view 
orientation for the hottest channel, 21lA. The temporally averaged maximum emission in 
this passband and orientation is 52 DNpix^^s^^. In these images we see that hot plasma 
fills most of the FC structure, without the small features that dominate the 171 A images. 



because the 21lA channel has a wide temperature response function (see, e.g., iLemen et al. 



201ll ). The hot material is evaporated from the chromospheric plasma by the footpoint 



heating. In both imag es, the core of the structure is bright, resemblin g the "chewy nougats 



observed by Yohkoh (IHudson et al.l Il999l : iHudson &: SchwennI |2000| ) . In contrast with the 



17lA images, however a diffuse halo surrounds the prominence in 21lA emission reaching a 



height of 40 Mm, resembling the "hot prominence shrouds" reported by lHabbal et al.l (120101 ). 
In the side views (Figs. [TTfe and [TTf ) the emission is more uniform, with no distinguishable 
bright structure in the core. In this orientation the time-averaged maximum emission is 
28 DNpix~^s~^. This demonstrates that the intensity distribution depends on the relative 
orientation of the FC to the LOS. In particular, the bright core appears when only the PIL 
is nearly aligned with the LOS. In the end view, the flux tubes in the core are closely aligned 
with the LOS, producing the strong core emission. In the side views (Figs. [TTb and [TTIF) . the 
21lA emission is more uniform and the core is darker because the cool plasma of the promi 



nence emits little in this bandpass, an effect also known as emissivity blocking (lAnzer et al. 



20071 ). In addition, the emission in this channel is more or less constant throughout the 
simulation time (see the accompanying movies) in both orientations. This indicates less 
dynamism of the hotter plasma than the cooler plasma. 



4.1. Differential emission measure 



In Figure [121 the differential emission measure (DEM) of our prominence model and 
overlying arcade is shown, taking into account all of the flux tubes of the filament channel 
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structure. This is averaged in space and time for each flux tube, and then combined over 
the whole volume. To avoid the contribution of the initial transient phase, we averaged 
the DEM from 5 hrs after t he onset of the localized he ating to the end of the simulation. 
The DEM is computed as in lKarpen fc Antiochod (120081 ). The result of this computation is 
normalized by a prescribed column depth, which is chosen to ensu re a good fit with the DEM 
values derived from prominence observations by SOHO SUMER (jParenti fc Vialll2007l ). The 
resulting column depth is ~ 10 km, two orders of magnitude below our previous results 
{h ^ 1000 km) because the present work includes approximately 100 times as many tubes. 
As we see in the figure, our model reproduces very well the shape of the observed quiet- 
prominence DEM between log T = 4.6 and 5.7. In addition, the location of the relative 
mi nimum of our profile is log T = 5.1, consistent with the observational value. As argued 



m 



Karpen &: AntiochosI (|2008[ ). the additional emission at higher temperatures (log T ^ 6) 



can be attributed to the unsubtracted back ground and foreground coronal contribution to 
the observed DEM of iParenti fc Viall (120071 ). Moreover, in our set of selected flux tubes we 
have considered only a small fraction of the overlying arcade loops that exist in a realistic 
structure. These loops of the overlying arcade contribute mainly to the warm and hot part 
of the DEM function because the condensations formed in those loops are small. Therefore, 
increasing the number of overlying arcade loops should yield better agreement between the 
computed DEM and the observed curve at temperatures higher than log T = 5.7. We plan 
to extend our model to include more of the overlying arcade. 



5. Discussion and conclusions 



We have developed a comprehensive three-dimensional model of the plasma dynamics 
and energetics of a prominence and its overlying arcade. We have found that condensations 
are formed in all of the flux tubes studied. There are two populations of condensations 
characterized by their cool mass. Repeated cycles of condensations forming and falling 
to the chromosphere produce condensations of masses below 7 x 10~^ g cm~^, and small 
sizes. Condensations that linger in the dipped portion of the tube have a constant growth 
rate, and reach large masses above this limit and large lengths. We have called the small 
condensations "blobs" and the large ones "threads". Blobs have a mean length of about 
1 Mm, while threads are larger condensations with a mean length of 37 Mm. The dominant 
contribution to the prominence mass comes from the threads. This mass grows at a constant 
rate of 39 x 10^^ g cm~^ hr~^. Assuming that the threads have a width of 200 km, the total 
mass of our model prominence is 1 x 10^^ g at the end of the simulation, in agreement with 
observed values for small prominences. When threads are formed, they oscillate around an 
equilibrium position; this motion is mainly horizontal (parallel to the solar surface), with a 



o 




4.5 



5.0 



5.5 



6.0 



6.5 



log(T) (K) 



Fig. 12. — Time- averaged differential emission measure vs. temperature for the whole system 
of flux tubes of the filament chan nel and the correspon ding overlying arcade (solid line). 
The observed best fit obtained by iParenti &: Viall (120071 ) for prominence section A_l (data 
courtesy of Parenti 2007 [private communication]) also is shown (dashed line). 



mean velocity for the entire set of threads of 34 km s~^ and a maximum value of 83 km s~^. 
Blobs move in all directions from horizontal to vertical, with a mean velocity of 51 km s~^ 
and a maximum of 92 km s~^. 

With our three-dimensional model we have constructed synthetic images including the 
orientation of the system and LOS integration effects. The synthetic Ha images reveal two 
trapezoidal-shaped prominence segments, each centered at a magnetic dipole and made up 
primarily of threads. These images clearly show the threads growing and the counterstream- 
ing motion occurring on adjacent threads. We have found that blobs are ubiquitous features 
in these images, moving along the field lines and falling to the chromosphere. Most of the 
blobs decelerate before entering the chromosphere. Blobs of the overlying arcade are the 
coronal rain, with primarily downward motions. Blobs are also produced in the sheared 
flux tubes associated with the prominence, but their motions are more horizontal than the 
coronal rain, consistent with observed counterstreaming. We have also computed synthetic 
images in two EUV channels, 17lA and 21lA. Most of the predicted 17lA emission comes 
from the PCTR - the warm plasma between the threads or blobs and the adjacent hot corona. 
The PC core is brightest in the end view because the concentration of threads and blobs is 
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highest there. In the zone immediately surrounding the core, only blobs are formed so the 
171A emission is collectively dimmer and more variable. In addition, long, bright, transient 
features occupy most of a tube just before a blob is produced, as the plasma cools through 
the 17lA passband. In 21lA, strong emission comes from the FC core in the end-view orien- 
tation, with a diffuse halo surrounding the prominence. This indicates that the cool material 
of the p rominence is surrou nded by hot material resembling the "chewy nougats" seen by 
Yohkoh (iHudson et al.lll999l) and the observed "hot prominence shroud" from recent eclipses 
( jHabbal et al.ll2010l ). We have computed the D EM and found very go od agreement with the 
observational DEM of a quiescent prominence (jParenti fc Viall 120071 ). The computed DEM 
curve reproduces the observed shape between temperatures log T = 4.6 and 5.7, and the 
observed position of the minimum at log T = 5.1. 

The ability of our model to reproduce many observed prominence characteristics lends 
strong support to the thermal nonequilibrium mechanism for producing the mass of the 
prominence, as well as the plasma dynamics. The agreement between model and observations 
also suggests that the fundamental structure of the filament channel is well described by the 
sheared arcade. Currently, we are extending the realism of our model in three directions, to 
allow quantitative comparison with prominence and cavity observations. First, we are greatly 
increasing the number of tubes that sample both the FC and the overlying arcade. Second, 
in this work have assumed that the footpoint heating is steady in time and is the same on 
all flux tubes. In an upcoming study we will investigate the effects of a heating function 
th at depends on the key physical parameters, according to the scaling laws summarized 



m 



Mandrini et al. (I2OO0I), as wel l as allow for temporally varying heating, consistent with 



nanoflaring (e.g.. iKlimchukI 120061 ). Finally, the emission in EUV coronal lines is significantl y 
affected by Lyman absorption by neutral H and He and by He^ (e.g., iKucera et al.lll998l ): 
these important effects will be included in future visualizations of our simulations. The large 
amplitude oscillations of the condensations in our prominence model will be addressed in an 
upcoming publication. 
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